Broken translational and time-reversal symmetry in superconducting films 
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We demonstrate that films of unconventional (d-wave) superconductor at temperatures T < 
0.43 T c can exhibit unusual superconducting phases. The new ground states beside the broken 
gauge and the point group symmetries can spontaneously break (i) continuous translational sym- 
metry and form periodic order parameter structures in the plane of the film, or (ii) time-reversal 
symmetry and develop supercurrent flowing along the film. These states are result of the strong 
transverse inhomogeniety present in films with thickness of several coherence lengths. We show a 
natural similarity between formation of these states and the Fulde-Ferrell-Larkin-Ovchinnikov state. 

PACS numbers: 74.81.-g,74.78.-w,74.20.Rp 
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Introduction. Properties of superconducting state in 
confined geometry will have defining influence on mini- 
mization of superconducting devices. This is especially 
true of an unconventinal superfluid with order param- 
eter (OP) that breaks more than one of the normal 
state symmetries 0, 0, HJ and which can be suppressed 
by interfaces, forming new non-uniform ground states. Q 
Such states can exhibit new broken symmetries, and 
currently are subject of broad theoretical and experi- 
mental investigation. Q For example, they can arise in 
bulk superconductors under influence of external field 
or pressure. For example, magnetic field breaks the 
time-reversal symmetry by inducing supercurrents, or, 
as in case of Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 
state0, 0], results in an oscillating OP that breaks con- 
tinuous spatial translation symmetry. In unconventional 
superconductors new ground states have been predicted 
to appear spontaneously even in the absence of exter- 
nally applied magnetic field. One origin of these states is 
the interaction of self-induced magnetic field, caused by 
distortion of the OP, with the OP itself. ||, @, 0, ED, EL1| 
Other current-carrying states arise near surfaces due to 
subdominant d + is pairing that breaks time-reversal 



symmetry. [13l . |14| Finally, new states with triplet struc- 
ture may appear in quantum wires. (l5l| Ultimately, the 
origin of these states lies in the app earance of Andreev 
states bound to interfaces, p, 11, 16[ 



In this Letter we propose existence of new ground 
states with non-trivially broken symmetries that require 
neither the self-induced magnetic field nor complex sub- 
dominant pairing. These states are induced by strong 
distortion of the OP shape. We consider this problem 
in the context of a singlet d-wave superconductor, when 
distortion of the OP is caused by confinement to a film 
or wire. We assume that the thickness of the film is 
ajustable, and that by varying it we are able to drive 
transitions between different ground states. The sim- 
plicity of this model makes the problem more transpar- 
ent, yet allows to draw general conclusions about connec- 
tion of inhomogeniety and spontaneous symmetry break- 
ing in unconventional pairing states, including superfluid 



3 He.[I3 

Free energy in a film. We consider a film of a d-wave 
superconductor with pairbreaking surfaces, such that the 
OP in the film is very non-uniform. This is achieved by 
orienting a gap node along the film (sketch in Fig. [I}. 
We assume a cylindrical Fermi surface, and OP solutions 
uniform along its axis (z-direction). Below the super- 
conducting bulk temperature T c , the OP can be contin- 
uously suppressed to zero by reducing the thickness D 
of the film, eventually reaching the normal state below a 
critical value of D* (T) . [lj| Our goal is to look for states 
that have more features than just the trivial suppression 
of the OP near the film's surfaces. 

We start with derivation of the Ginzburg-Landau 
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FIG. 1: (Color online) A d-wave superconductor (SC) con- 
fined to a film with pairbreaking surfaces. (a) The SC- 
normal(N) transition as a function of temperature for sev- 
eral (fc-modulated states. At low temperature, transition into 
uniform (q x = 0) SC state shows re-entrant behavior and is 
preempted by transition into a state with finite ^-modulation; 
(b) the envelope of various q x curves gives the maximal con- 
finement vector Ql(T,Q x ) = n/D*, plotted together with 
the corresponding modulation vector Q X (T). The coherence 
length £ = hvf/2-KTc. 
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free energy that describes this transition. For this we 
solve the microscopic quasiclassical equations [IH [2(J for 
the diagonal, g, and off-diagonal (anomalous) compo- 
nents, / and /', of the quasiclassical Green's function, 
<?(R, p; e, n ). Complete set of equations also includes the 
normalization condition and symmetry relations, 



^v(p). V R ]/ = i<? A(R,p) ; 



(1) 



9 2 -ff 



/'(R, p; e m ) = /(R, p; -e m )* 



Here v(p) = Vff> is the Fermi velocity at point p on the 
Fermi surface and e m = 7rT(2m + 1) - Matsubara en- 



J 



ergy. We assume a separable OP, and expand the spatial 
part in plane waves, A(R, p) = y(p) £ q A q e iqR . The 

normalized basis function is y(p) = \/2 s m20p, and the 
plane wave amplitudes satisfy self-consistency equation, 



y(P) /q(p;6m) 



^A q y(p) 



(2) 

that follows from minimization of the free energy func- 
tional, SaF/SA^ = 0, and brackets denote angle average, 

(•••)= / "2^r • • ■ ■ Solving Eqs.([r| for / up to third order 
in A q , we use self-consistency ©, to find the free energy, 



AF = Y, /(T,q)|A q | 2 + | A'(T, qi ,q 2 ,q 3 ,q4) A qi A q2 A q3 A c 



qi+q2=q3+q4 



/,r, q ,.i„l-2,r^ R o^,p)(_L--- 



where ?7 q = -v • q , 



tf(T,q ljq2 ,q 3 ,q4) = 27rT ^ lR e /y\p) 



e m >0 



£rn + ^fa + Vq 2 + Vq a + T?q 4 )/4 

(e m + «T? qi )(£m + i?7q 2 )( e m + ir)q 3 ){e m + z?y q4 ) 



(3a) 
(3b) 
(3c) 



This functional describes second-order transition from 
uniform normal state into general, modulated supercon- 
ducting state. For analytic analysis and to make use of 
the fewest number of plane waves in the OP we consider 
perfectly specular film surfaces with boundary condition 
f(x,±D/2,p]E m ) = f(x, ±D/2,p;e m ), connecting the 
incoming, p, and mirror-reflected, p = p 2y(yp), 
trajectories at surfaces y = ±D/2. From the self- 
consistency it follows that the OP at the surfaces is zero 
due to d-wave symmetry. This ensures the OP form 
A(R) = J2 Qao A gx e l<?x:r cos ny/D, and thus fixes the y- 
wave number q v = Q y = n/D. Note, that fixed Q y in 
? /q = ( v xQx + v yQy)/2 (Eqs.[3]) plays the same role as the 
external magnetic field in ?y q B = v • q/2 + [iB that enters 
similar equations describing instability into FFLO phase 
in Pauli-limitcd superconductors. 21 

Phase diagram of superconducting state in a film. We 
first find the boundary of the superconducting phase, and 
thus look for the largest Q* y {T) (narrowest strip D*(T)), 
where superconductivity first appears. This instability 
is given by condition I(T, q x , Q*) = 0, and its solution, 



Q*(T,q x ), is shown in Fig.[T^, for various q x states. The 
q x = state, uniform along the film, at low temperature 
is preempted by states with finite modulation q x . For 
any given T we determine the optimal q x (denoted by 
Q X (T)) that gives the largest Q* (see Fig. [Tb); and find 
that x-modulated states are possible below T* ~ 0.43 T c . 

Next, we determine the relative stability of various 
non-uniform states near this transition. For each in- 
stability mode Q(T) = (Q x ,Q y ) there are three other 
degenerate modes obtained from this one by reflec- 
tion of x and/or y coordinates. We consider two 
principal states obtained by different combinations of 
these modes: (i) with two opposite g x -componcnts, cor- 
responding to the amplitude oscillations of the OP, 
which includes all four degenerate states, qi,2.3,4 = 

{(Qx, Qy) , (Qx, ~Qy) , (~Qx, -Qy) , (-Qx, Qy)}', and 

(ii) with one g x -component, that gives a modulation of 
the OP phase along the film and therefore a superflow, 
Ql,2 = {(Qx, Qy) , (Qx, -Q y )}- 

Near the second-order transition. 



AF[Q X ,-Q X ] 



2/ 2 (T,Q) 



K 1 + 2K 12 + 2K 13 + 2K U + 2K 123 4 
2 / 2 (T, Q) 



AF[Q X ] = - 



2(Ki + 2K 12 ) 



for A(R) = Ax(e iQa:X + e ~ lQ * x ) cosQ y y , (4a) 
for A(R) = A 2 e iQ ^ cos Q y y . (4b) 
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FIG. 2: (Color online) Free energy density of a d-wave su- 
perconductor confined to a film, as a function of the inverse 
thickness for T/T c = 0.1. In films with D < D* = tt/Q* the 
states with superflow have lower energy than the q x — state. 
The inset shows self-consistently calculated order parameter 
profile, A(y), and the current density, J s ,x(y), for small q x 
and D < D*. The current has anomalous (paramagnetic) 
contributions near the edges due to Andreev bound states, 
and Qy marks the point where the average superfluid density 
in the film, p s =. lim J® x /q x , vanishes. 



Here we use notation, K\ = K(qx, qi, qi, qi), if 1234 = 
^"(qii q2,q3,q4), an d for pairs of wave vectors, = 
qj> qii q? )- At low temperatures we calculate these 
coefficients analytically and find that -K13 term, corre- 
sponding to the two opposite wave vectors ±Q (Q \ 
node), diverges as 1/T. This makes the current-carrying 
state with spontaneously broken time-reversal symme- 
try be the lowest in energy. 22| Numerical evaluation of 
all if's shows that the state with a current has low- 
est energy at all T < T*. We also find that transi- 
tion from normal to superconducting state is always sec- 
ond order in the film, to be contrasted with bulk Pauli- 
limited superconductors, where A 4 -term coefficient be- 
comes negative and transition becomes first order at low 
temperatures. [IE 0| In a film, even for q x = 0, we have 
a wavevector pair ±Q y 7^ and K12 > 0, which guaran- 
tees positive sign of K\ + 2K\% in Eq. (|4b|) . even though 
K\ < at low temperatures. 

To fully describe the phase diagram and elucidate the 
structure of the new phases, we determine the transi- 
tion between the new phases and the q x = 'uniform' 
condensate deep inside the SC state. We start with the 
state that breaks translational symmetry and forms pe- 
riodic modulations of |A(R)| in the film's plane. In this 
case the general form of the OP is similar to one in the 
FFLO problem. 25, H(| Near the normal state instabil- 
ity A(x, y) ~ cosQ x x cosQ y y, but becomes less har- 
monic and more domain- like structured along x, as we 
increase film thickness (Q y decreases), The distance be- 
tween neighboring domains grows and at some critical 
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FIG. 3: (Color online) The phase diagram of a d-wave su- 
perconductor in a film. The shaded region marks the ground 
state with flowing current. Suppression of this state would 
open access to the state with broken translational symmetry. 



ate states at x = ±00 with opposite OP profiles ±A(y), 
disappears (square-symbol line in Fig. [3]). Note, that this 
transition occurs below unphysical q x = line (thin dots 
in Fig. [3]), as compared with the FFLO problem 25|, 2(|: 
there is no first order transition above q x = line in this 
case - thus the inhomogeneous phase must cover the en- 
tire thermodynamically unstable region. Energetically, 
the state with the modulated amplitude of the OP gains 
energy compared with the energy of 'uniform' state from 
the reduction of pairbreaking at the domain walls and 
the redistribution of the Andreev bound states. [l7| 

Next, we consider state with phase modulation (su- 
perflow), A(R) = A cos Q y y e z< ^ lX , that carries super- 
current, J s ,x{y) = 2N f TJ2 em (v x {p)g{y,p;s m )). We 
write down a more general form of the OP, A(x,y) = 



A 9x (y) exp(iq x x) , and self-consistently determine its am- 



plitude profile and the associated free energy density [2 

fD/2 
-D/2 



n D (q x )= / ^AF D (q x ,y), 

J-D/2 U 



(5) 



Q y the last domain wall, that separates two degener- 



as a function of q x and the film thickness, D. Key fea- 
tures of this calculation are shown in Fig. O In the main 
panel we plot QP as a function of inverse film thickness 
for several q x . States with finite superflow are stabilized 
in films thinner than D*(T) = n/Q*(T). Moreover, their 
stability region extends beyond that of q x = state. 
Since finite q x induces supcrcurrent, it can be the ground 
state only when the total current in the film disappears 
= P? (<&)& = dn D {q x )/dq x = 0.EH Upper critical 
width D* is determined by vanishing average superfluid 
density, (0) = 0. These conditions are possible to sat- 
isfy due to backffowing anomalous surface currents car- 
ried by Andreev bound states (inset of Fig. [2]). HI, EJ, 0] 
The complete phase diagram is presented in Fig. [3] 
Below T* ~ 0.43 T c , two new ground states are possible 
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k + - k = 2 |iB / v ( 



FIG. 4: (Color online) (a) inhomogeneous state of FFLO type, 
when the pair-forming states with opposite spins belong to 
different Fermi surfaces, split by magnetic field; (b) state with 
spontaneous current along x. Fermi surface is shifted by fixed 
center-of-mass flow ±Q y . Paired states can have additional 
i-momentum Q x to minimize pairbreaking. 

in superconducting films. The state with spontaneous 
current and broken time-reversal symmetry takes a large 
part of the phase space. Under considered conditions 
the state with modulated amplitude of the OP lies inside 
the stability region of the current-carrying state and is 
not realized. However the relative energies of the two 
states may be affected, e.g. by surface roughness. Also 
note, within numerical precision the lower instability line, 
D*(T), is nearly straight and extrapolates to the origin. 

Qualitative picture for appearance of the longitudinal 
modulations of the OP in the film is illustrated by anal- 
ogy with the FFLO state in Fig. g] In the FFLO state 
modulations of the OP arise to minimize pairbreaking 
caused by pairing of electrons across Zceman-split Fermi 
surfaces, ± [iB, Fig. 0Ja). In the strongly inhomo- 
geneous state in the film, the role of magnetic field is 
played by the fixed confining wave vector Q y that pro- 
duces superflow across the film and shifts the Fermi sur- 
face vertically, £k ± v y Q y /2, Fig. 0Jb). Pairs are formed 
with additional shift Q x (so that Q l/f node) to minimize 
pairbreaking by the superflow, c.f. Ref. I29I . 

Finally, we remark on the observability of these 
states. In superconductors diamagnetic coupling to (self- 
induced) magnetic field will modify the phase digram, 
Fig. [3J Previous studies of semi-infinite (1/D — > 0) 
systemlUGilEl 

show that spontaneous surface current 
appears already at finite temperature and this indicates 
that the phase space for q x ^ state may be enlarged. 
Also, modulated states in confined geometry are less sen- 
sitive to the disorder and would persist much longer in 
dirty samples compared with surface states. [13] In fact, 
we find that complete suppression of T* — > by impuri- 
ties happens only when the corresponding suppression of 
T c is 60% (mean free path I ~ 5£o)- It is also reasonable 
to expect that as long as the OP has significant gradient 
across the film, the new superconducting phases are only 
slightly affected by surface roughness. 

Conclusions. We have studied behavior of a d-wave 
superconductor in a film geometry. We find that large 



gradients of the order parameter across the film, can be 
- quite countcrintuitively - 'relieved' by producing addi- 
tional modulation along the film. This means that in 
confined geometry this pairing system may undergo a 
transition into a state with broken translational or even 
time-reversal symmetry. Similar behavior occurs in a 
completely different pairing state (triplet, p-wave) in 3 He 
films [T3l and we suggest that these states are common to 
unconventional pairing systems, although their exact na- 
ture may depend on the structure of the OP. The new 
ground states are more robust in film or wire geometry 
than in bulk or semi-infinite systems and should be ob- 
servable in superconductors and neutral superfluid 3 He 
in confinement. 
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